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Resumen 


La presión real de vapor de agua es una variable básica para estimar la 
evapotranspiración de los cultivos, uno de los componentes del ciclo hidrológico; 
sin embargo es difícil y cara de medir de forma directa, por lo que se recurre en 
la práctica a estimaciones basadas en la temperatura y relaciones sicrométricas. El 
objetivo del presente trabajo fue realizar una comparación de diferentes métodos 
convencionales para el cálculo de la presión real de vapor y compararlos con las 
estimaciones realizadas con dos tipos de redes neuronales artificiales: feedforward 
backpropagation y radial basis function. Se usaron datos meteorológicos de cuatro 
estaciones del Distrito 075, localizadas en el Valle del Fuerte, al norte del estado 
de Sinaloa, México. Los resultados indican que la red neuronal artificial tipo radial 
basis function (escenario E4) mostró ser el mejor método en la estimación de la 
presión actual de vapor de agua. 


Palabras clave: humedad atmosférica, déficit de presión de vapor, redes neuronales 


artificiales. 


Introducción 


La presión real de vapor de agua es un elemento 
importante en la existencia de la humedad 
circundante en la atmósfera. Determinar el 
déficit de presión de vapor resulta necesario 
para el cálculo de la evapotranspiración (Allen 
et al., 1998), variable relacionada con los re- 
querimientos de agua de los cultivos (Jensen et 
al., 1990), y que influye finalmente en la lámina 
de riego durante la programación de riegos de 
los cultivos; de igual manera, la presión real de 
vapor de agua representa un elemento básico 
como variable en la modelación de crecimiento 
de cultivos y, como lo indica Baker (2005), no 
existen muchos métodos directos para medir la 


presión o densidad de vapor. La presión de vapor 
sólo puede ser medida de forma directa por 
medios ópticos, que son caros para la mayoría 
de las aplicaciones agrícolas; en consecuencia, 
su estimación práctica se realiza a partir de la 
temperatura y de relaciones sicrométricas. Se 
han propuesto diversos modelos para estimar 
la presión real de vapor, como los expuestos 
por Allen et al. (1998), que plantean el uso de la 
temperatura de punto de rocío; para su cálculo 
también se han propuesto diversos métodos, 
como la ecuación de Murray (Howell y Dusek, 
1995) o la que proponen Linsley et al. (1998). 
Para la propia presión de vapor actual, Allen et 
al. (1998) presentan una ecuación que requiere 
la temperatura mínima y la humedad relativa 
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máxima. Sadler y Evans (1989) demostraron 
con 15 métodos diferentes de cálculo del 
déficit de presión de vapor, que el error varió 
de un 80% en subestimación al 100% en la 
sobreestimación de la evapotranspiración, 
comparados con el mejor método para estimar 
el déficit de presión de vapor de agua. 

Las redes neuronales artificiales (RNA) 
son un modelo que aprende por medio del 
entrenamiento y captura relaciones funciona- 
les entre los datos, incluso si las relaciones 
subyacentes no son conocidas o difíciles de 
describir, y es en este sentido que las RNA 
pueden ser tratadas como un método estadístico 
multivariado no lineal y no paramétrico (Zhang 
et al., 1998). Estas características representan 
ventajas que colocan a las RNA como un 
excelente método para predecir y/o estimar 
variables meteorológicas. 

Los autores no encontraron investigacio- 
nes para estimar la presión de vapor real de 
agua usando redes neuronales artificiales 
(RNA), pero se han realizado diversos trabajos 
para estimar el déficit de presión de vapor, la 
humedad relativa y la temperatura de punto 
de rocío con RNA (Mittal y Zhang, 2003; Shank 
et al., 2008a,b). Bialobrzewski (2008) estima la 
humedad relativa en el aire entrenando una 
RNA feedforward multicapa con retraso en el 
tiempo, con el algoritmo backpropagation para 
estimar 16 datos en intervalos de tres horas, es 
decir, una predicción de la humedad relativa 
en las próximas 48 horas. He et al. (2007), por su 
lado, predicen por medio de una red neuronal 
la humedad relativa y la temperatura del 
aire en un invernadero. Shank et al. (2008a,b) 
estimaron la temperatura de punto de rocío, 
paso previo en el cálculo de la presión real de 
vapor de agua, por medio de una RNA tipo 
backpropagation de forma horaria. 

El objetivo de la presente investigación 
fue llevar a cabo una comparación de la pre- 
sión real de vapor de agua estimada con los 
métodos convencionales de Linsley, Murray y 
Bosen contra las estimaciones realizadas con 
un modelo de red neuronal artificial multicapa 
tipo feedforward backpropagation y otro radial 
basis function, y la validación de estos modelos. 


Marco teórico de la RNA multicapa 
feedforward backpropagation 


Las redes neuronales artificiales, de manera 
reciente han probado ser excelentes herra- 
mientas para aproximar sistemas no lineales. 
El algoritmo usado en el modelo feedforward 
multicapa de red neuronal artificial fue el 
backpropagation. Este algoritmo empieza con 
la inicialización, la presentación de ejemplos 
para el entrenamiento y para cada ejemplo se 
ejecutan los dos pasos siguientes de acuerdo 
con Haykin (1994): 


Paso 1. El cálculo hacia delante. Sea un ejemplo 
de entrenamiento en la época denotado por 
[x(m), d(n)], con el vector de entradas x(n) 
aplicado a los nodos de la capa de entradas y el 
vector de respuestas deseado d(n) presentado al 
nodo de la capa de salida. El nivel de actividad 
interno de la red 0/An) para la neurona ¡en la 
capa l es: 


Pp 
9 (1) = ws (my) E 
i=0 


donde y/*(n) es la función señal de la neurona 
¡ en la capa previa ! — 1 en la iteración 1; 10,M(n) 
es el peso sináptico de la neurona ¡ en la capa 
[, que es alimentada de la neurona ¡en la capa 
1 — 1. Para i =0, se tiene y" (1) =-1 y yo 
(n) = as, donde 9/1) es un umbral aplicado 
para la neurona j en la capa !. Asumiendo el 
uso de la función logística para la no linealidad 
sigmoidea, la función señal (de salida) de la 
neurona ¡ en la capa l es: 


Dn - 1 
UN (n) 1+expoXm)) (2) 


Si la neurona j está en la primera capa 
escondida (esto es, | = 1), se tiene: 


ym) =x,(1) (3) 
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donde x;(1) es el j-ésimo elemento del vector de 
entrada x(n ). Si la neurona j está en la capa de 
salida (1 = L), se tiene: 


ym) =0/(1) (4) 
de aquí, se calcula el error de la señal: 
e/(1)=d/(1)-0,(1) 65) 


donde d¡(n) es el j-ésimo elemento del vector de 
respuestas deseado d(n). 


Calcula los 
Ys (gradientes locales) de la red por el 
procedimiento hacia atrás, capa por capa: 


Paso 2. Cálculo hacia atrás. 


98m) = m0, (1)[1-o/(m)] (6) 


para la neurona ¡ en la salida de la capa L: 


9 (n)= ym] - IDO (7) 
k 


para la neurona j en la salida de la capa l. 

De aquí se ajustan los pesos sinápticos de 
la neurona en la capa l de acuerdo con la regla 
delta generalizada: 


(1 ( 
wi (+1) = w/, (11) 


+ asp (m) se wn S D]+ my (m) (8) 


donde mn es el parámetro tasa de aprendizaje 
y a. es la constante de momento; y el último 
componente del algoritmo backpropagation es la 
iteración donde se presenta nuevas épocas de 
ejemplos de entrenamiento a la red hasta que 
los parámetros libres de la red estabilicen sus 
valores y el error medio cuadrado (E, ) calculado 
sobre el conjunto entero de entrenamiento 
tenga un valor mínimo o pequeño aceptable 
(Haykin, 1994). 


Marco teórico de la red neuronal artificial 
multicapa radial basis function (RBF) 


Una alternativa viable a la alta no linealidad 
en los parámetros de la red neuronal es la 
red radial basis function (RBF), que puede ser 
considerada como un caso especial de una red 
de dos capas, la cual es lineal en los parámetros 
debido al fijamiento de todos los centros RBF y 
no linealidades en la capa oculta. De esta forma, 
la capa oculta ejecuta una transformación no 
lineal fija con parámetros no ajustables y mapeo 
al espacio de entradas en un nuevo espacio. 
La capa de salida entonces implementa una 
combinación lineal sobre este nuevo espacio y 
los únicos parámetros ajustables son los pesos 
de esta combinación lineal. Estos parámetros 
pueden, por lo tanto, ser determinados usando 
el método lineal de mínimos cuadrados. En la 
práctica, los centros de la RNA RBE son elegidos 
con frecuencia como un subconjunto de los 
datos. La mayoría de los trabajos publicados 
simplemente asumen que los centros se 
seleccionan de manera arbitraria a partir de los 
datos (Chen et al., 1991). 


Materiales y métodos 
Área de estudio y obtención de datos 


Se utilizaron datos meteorológicos de hume- 
dad relativa (%), temperaturas máximas y 
mínimas (*C), obtenidos de la base datos de la 
red agroclimática automatizada del Valle del 
Fuerte, localizada en el Distrito de Riego 075, 
ubicada en la vecindad de la ciudad de Los 
Mochis, Sinaloa. Tales datos fueron de cuatro 
estaciones meteorológicas (figura 1): 3843 II-2 
Ruiz Cortines, 3546 II-3 Batequis, 3765 I11-1 AC 
Santa Rosa 1 y 9610 III-1 AC Santa Rosa 2, cuyas 
coordenadas geográficas de latitud, longitud y 
altitud (msnm) son, respectivamente: 25* 39' 
157, 108” 45' 20”, 31; 25* 45' 49”, 1087 48' 41”, 32; 
25” 45' 03”, 108* 57' 21”, 40; y 25* 51 16”, 108* 
52' 03”, 61. Estos datos se preprocesaron todos 
los días. Comprendieron del periodo de abril 
de 1997 a mayo de 2001, generando vectores 
de 1 484 datos diarios, los cuales se utilizaron 
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para el ajuste de los modelos convencionales, y 
para el entrenamiento, validación y prueba de 
las RNA. Otro conjunto de datos del periodo de 
junio a diciembre de 2001 (229 días) se utilizó 
para realizar, con los diferentes modelos, una 
predicción (validación) de la presión real de 
vapor de agua. Los datos se preprocesaron en 
hojas de cálculo de Excel de Microsoft Office 
2007 y el software usado para el entrenamiento 
de las RNA fue el Matlab 7.0 (Demuth et al., 
2008), con sus utilerías (toolbox) para redes 
neuronales. 


Implementación de la redes neuronales 
artificiales feedforward backpropagation y 
radial basis function 


En el entrenamiento de la RNA feedforward 
backpropagation se consideraron diferentes con- 
juntos de variables de entrada con una sola 


capa oculta y una sola neurona en la capa de 
salida; se utilizó el algoritmo de entrenamiento 
Levenberg-Marquardt, y la función de trans- 
ferencia tangente-sigmoidea (tansig) en la 
capa oculta y en la capa de salida; todos los 
entrenamientos convergieron en menos de 
cien épocas. Los datos usados (1 484) se ocu- 
paron de la siguiente manera: la mitad para 
el entrenamiento, una cuarta parte para la 
validación y la otra cuarta parte para la prueba 
de la red neuronal artificial. Como parte de la 
propia arquitectura de la red, se consideraron 
combinaciones de las variables de entrada, de 
la temperatura máxima, temperatura mínima, 
humedad relativa máxima y humedad relativa 
mínima, con cuatro neuronas en la capa oculta 
y como salida, la presión real o actual de 
vapor de agua; previo a los entrenamientos se 
escalaron las entradas y salida a un rango de 
[-1, 1]. 


í Santa Rosa 2 


ap 


a > 
Batequis 


Leyenda 


dis Estaciones 


[] Distrito de Riego 075 


Figura 1. Región de estudio. 
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En la RNA radial basis function, el extendido 
o spread elegido fue de cien en todos los 
entrenamientos; el código en matlab para esta 
función es “net = newrbe (P, T, spread)” donde se 
crean tantas neuronas radbas como vectores de 
entradas existen, y donde P son los vectores de 
entrada y T es el target o vector de salida. 


Modelos convencionales 
Uno de los modelos utilizados para estimar 


la temperatura del punto de rocío es con la 
expresión de Murray (1967), expresada como: 


es 237.3 
. -1 
HR 
In| — 
(700) E E 
17.27 | |237.3+T (9) 
Donde T,,.., es la temperatura del punto 


de rocío (*C); T, temperatura media (*C); HR, 
humedad relativa media (%). 

Otra forma alterna de estimar la tempera- 
tura del punto de rocío está dada por la 
expresión mencionada por Linsley et al. (1998), 
y es: 


T 


rocío — 


T-(14.55+0.114*T)*X 


-[(2.5+0.007*T)* X]'-(15.9+0.117*T)*X"* (10) 
Donde: 
Otra ecuación para estimar la temperatura 


de punto de rocío está dada por la expresión de 
Bosen (1958): 


1/8 
Caco" (595) *[112+0.9*T]+(0.1*T)-112 (11) 


el valor T 


rocío” 


calculado con las expresiones (9), 
(10) y (11), se sustituye en la ecuación siguiente 
para estimar la presión real de vapor le kPa) 
(Allen et al., 1998): 


17.27T 


rocío 


e,= e(T a LAR 
(+27) 


rocío 


) =0.611exp 


| (12) 


Allen et al. (1998) indican que cuando 
no se tiene el valor de T 
la T 


mín 


se puede utilizar 
(*C), ya que la relación T 


rocío 


rocío 
- se 
min 
mantiene para localidades en donde un cultivo 
bien regado cubre el suelo de la estación. Sin 
embargo, como mencionan Allen et al. (1998), en 
regiones áridas, donde el aire no está saturado 
cuando se presenta la temperatura mínima, es 


posible que T 


mín 


pueda ser mayor que T,_., y 
se requiera una calibración local para estimar 
la temperatura del punto de rocío. En estas 
situaciones, el valor de T,. 
(12) puede sustituirse por T 
restándole de 2 a 3 *C. Se realizaron las 


estimaciones de la presión real de vapor de 


en la ecuación 
O aproximarse 


mín 


agua (e,) de acuerdo con este procedimiento. 


Modelos utilizados en la estimación de la 
presión real de vapor de agua 


Se estimó la presión real de vapor por diferentes 
métodos, cada método fue un modelo. Para el 
modelo 1 se consideraron los entrenamientos 
por medio de la red neuronal artificial (RNA) 
multicapa feedforward backpropagation, y para 
el modelo 2, la RNA radial basis function con la 
metodología teórica propia indicada para redes 
neuronales artificiales (Park y Sandberg, 1991; 
Chen et al., 1991). 

Con la ecuación (12) se estimó la presión 
real de vapor para los modelos del 3 al 8 (con- 
vencionales), lo único que cambia es el valor 
que se sustituye en lugar de la temperatura 
rocío. Los modelos 3, 4 y 5 utilizan: T,. -2,T,. 
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-3 y Ty respectivamente. Los modelos 6, 7 y 8 


calculan primero la T_. con las ecuaciones de: 


Murray (ecuación (9)), Linsley (ecuación (10)) 
y Bosen (ecuación (11)), respectivamente y se 


utilizan en la ecuación (12). 
Datos observados 


Según Allen et al. (1998), la presión real de 
vapor (e,) se puede estimar de la humedad 
relativa con la ecuación siguiente: 


A AR máx +e (Tus) BR 


100 (13) 


donde e” (T,_..) y e” (T,.) son la presión de 
vapor a saturación (kPa), en función de la 
temperatura mínima (T,.) y temperatura 
máxima (T',._). La ecuación (13) es la expresión 
que de acuerdo con Allen et al. (1998) es la 
que mejor se aproxima al valor medido u 
observado de la presión real de vapor de agua; 
como se careció de datos directos de la presión 
real de vapor, se asumieron como medidos los 
calculados por la ecuación (13) y contra estos 
datos se compararon los datos de presión real 
de vapor de agua estimados por los demás 
modelos. 


Índices de evaluación estadística de ajuste de 
los modelos 


Se obtuvieron los siguientes índices estadísticos 
de prueba para los diferentes modelos: el 
error estándar promedio o raíz cuadrada del 
cuadrado medio del error (RMSE), calculándose 
con la siguiente ecuación: 


12 


1 
N 
Y (e, de y 
RMSE=| L2— (14) 


el error medio (MBE), llamado también sesgo o 
desviación, calculado con el fin de caracterizar 
la bondad de cada uno de los modelos, dado 
por la ecuación siguiente: 


N 
] (a; -+t;) 
MBE= Pili) (15) 
el coeficiente de correlación, R, dado por: 
N E e 
Bn DO - 1 Nt,-t) 
1/2 
N 0 N Luco 
(3:07) [6D ) (16) 


donde el coeficiente de determinación está 
dado por R?; a es el dato estimado por el 
modelo; t, el dato que se asume como medido; 
N, el número de observaciones o estimaciones; 
a, el promedio de los datos estimados por el 
modelo; £, el promedio de los datos del modelo 
que se asumen como medidos (Cai et al., 2007; 
Pereira, 2004). 


Resultados y discusión 


Al analizar las figuras 2, 3, 4 y 5, se observa 
que los datos estimados con el modelo 1, al 
compararlos con los datos observados, están 
alrededor de la recta 1:1 en las cuatro estaciones. 
El ajuste del modelo 2 presenta una relación 1:1 
en todas las estaciones. El modelo 3 (T,.. - 2) y 
el 4 (T,, - 3) subestiman los valores observados 
en Ruiz Cortines y Santa Rosa 1; en las otras 
dos estaciones, los puntos se distribuyen mejor 
alrededor de la recta a 45”, pero se comporta 
mejor el modelo 3. El Modelo 5 (T,. ) presenta 
una mejor dispersión de sus puntos en las 
estaciones Ruiz Cortines y Santa Rosa 1, y en 
las otras dos se observa una sobreestimación 
entre los valores de tres a cuatro (kPa). Los 
modelos 6, 7 y 8, que utilizan temperatura del 
punto de rocío, presentan una sobreestimación 
en todas las estaciones (está más acentuada en 
la estación Santa Rosa 2), aunque sus puntos 
son menos dispersos, comparados con los 
modelos 3, 4 y 5. 

A partir de este análisis previo se obtuvieron 
los estadísticos de prueba: raíz cuadrada del 
cuadrado medio del error (RMSE), el error 
medio sesgado (MBE), así como el coeficiente 
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Figura 2. Comparación entre los valores estimados y los observados (kPa) por los diferentes modelos 


de la presión real de vapor de agua para la estación Ruiz Cortines. 
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Figura 3. Comparación entre los valores estimados y los observados (kPa) por los diferentes modelos 


de la presión real de vapor de agua para la estación Batequis. 
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Figura 4. Comparación entre los valores estimados y los observados (kPa) por los diferentes modelos 


de la presión real de vapor de agua para la estación Santa Rosa 1. 
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Figura 5. Comparación entre los valores estimados y los observados (kPa) por los diferentes modelos 
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de determinación (R?), la pendiente (m) y el 
intercepto (b) de la relación entre los datos 
observados y los estimados de la presión real 
de vapor para las diferentes estaciones para los 
modelos 1, 2, 3, 5, 6, 7 y 8, y debido a que el 
desempeño del modelo 4 es mucho menor que 
los modelos 3 y 5 en cuanto a la estimación de 
la presión de vapor real de agua, se omiten los 
estadísticos de prueba para dicho modelo. 

En los cuadros 1 y 2 se presentan los 
parámetros (m y b) de los modelos 1 y 2, asícomo 
los estadísticos de su ajuste en cuatro escenarios 
o casos El, E2, E3 y E4, respectivamente; en 
ambos casos se observa que el ajuste de los 
cuatro escenarios es bueno y mejora conforme 
aumenta el número de las variables de entrada 
(la diferencia entre casos depende del número 
de variables meteorológicas usadas). En el 
cuadro 1, la nomenclatura del escenario E4 


Tu Toy HR, HR) (4 x 4 x 1) significa 
que se utilizaron las variables indicadas dentro 
de las llaves (temperatura mínima, máxima, 
humedad relativa mínima y máxima), que son 
cuatro variables de entradas, cuatro neuronas 
en la capa oculta y una salida, y es en este 
escenario E4, en donde se presenta un ajuste 
muy bueno, como se había observado en las 
figuras 2 a 5. 

Al comparar los estadísticos, R? y RMSE 
del cuadro 1 versus el cuadro 2, se tiene que: 
en el El, el modelo 1 es ligeramente superior 
en las tres primeras estaciones por milésimas. 
De igual manera, los escenarios E2 y E3 del 
modelo 1 son ligeramente mejores que los 
del modelo 2 y en el último escenario (E4), 
fue el modelo 2 el que presentó mejor ajuste 
estadístico. En el cuadro 3 se observa que los 
modelos 3 y 5 tienen similares valores de R?; los 


Cuadro 1. Estadísticos de prueba de ajuste entre los datos observados de presión real de vapor y los datos estimados 


por los diferentes entrenamientos de RNA feedforward backpropagation. 


Modelo 1. Arquitectura (Multilayer Perceptrón) m b R? RMSE MBE 
El (To...) 
Ruiz Cortines (1 x 4 x 1) 0.9455 0.1095 0.9409 0.1906 0.0022 
Batequis (1x 4x1) 0.9293 0.1324 0.9291 0.2081 0.0012 
Santa Rosa 1 (1x 4x1) 0.9303 0.1397 0.9318 0.2102 0.0009 
Santa Rosa 2 (1x4 x 1) 0.9226 0.1450 0.9193 0.2092 0.0005 
E 
Ruiz Cortines (2 x 4 x 1) 0.9534 0.0914 0.9487 0.1776 -0.0002 
Batequis (2 x 4 x 1) 0.9370 0.1217 0.9403 0.1912 0.0048 
Santa Rosa 1 (2x 4x1) 0.9366 0.1253 0.9391 0.1985 -0.0009 
Santa Rosa 2 (2x 4x1) 0.9356 0.1215 0.9364 0.1854 0.0011 
E3 (To... BR.) 
Ruiz Cortines (2 x 4 x 1) 0.9623 0.0717 0.9553 0.1659 -0.0025 
Batequis (2x 4x1) 0.9627 0.0713 0.9641 0.1481 0.0022 
Santa Rosa 1 (2x 4x1) 0.9514 0.0967 0.9533 0.1738 -1.85E-05 
Santa Rosa 2 (2x 4x1) 0.9612 0.0728 0.9616 0.1444 0.0004 
HA, BR. E) 
Ruiz Cortines (4 x 4 x 1) 0.9983 0.0032 0.9960 0.0152 -0.0002 
Batequis (4x 4x1) 0.9991 0.0015 0.9996 0.0142 -0.0002 
Santa Rosa 1 (4x 4x1) 0.9989 0.0018 0.9996 0.0157 -0.0004 
Santa Rosa 2 (4x 4x1) 0.9989 0.0021 0.9992 0.0206 -4.35E-05 


R?: coeficiente de determinación; RMSE: raíz cuadrada del cuadrado medio del error; MBE: error medio sesgado; m: pendiente; b: intercepto; 


T temperatura mínima; T”_. : temperatura máxima; HR 


mín máx" 


¡-ésimo. 


mín" 


: humedad relativa mínima; HR_. : humedad relativa máxima: E; escenario 


máx" 
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Cuadro 2. Estadísticos de prueba de ajuste entre los datos observados de presión real de vapor y los datos estimados 
por los diferentes entrenamientos de RNA radial basis function. 


Modelo 2. RNA Radial basis function m b R? RMSE MBE 
El (To...) 
Ruiz Cortines 0.9393 0.1194 0.9393 0.1931 6.47E-06 
Batequis 0.9254 0.1384 0.9254 0.2136 2.76E-06 
Santa Rosa 1 0.9302 0.1391 0.9301 0.2127 7.08E-06 
Santa Rosa 2 0.9224 0.1450 0.9224 0.2050 -1.3E-06 
EPIA: 
Ruiz Cortines 0.9461 0.1061 0.9461 0.1819 -1.27E-07 
Batequis 0.9344 OBT2TE 0.9343 0.2002 5.77E-07 
Santa Rosa 1 0.9389 0,1217 0.9389 0.1989 -2.94E-07 
Santa Rosa 2 0.9282 0.1342 0.9281 0.1972 1.18E-06 
ES Tas HR.) 
Ruiz Cortines 0.9538 0.0907 0.9539 0.1683 -4.61E-06 
Batequis 0.9633 0.0680 0.9633 0.1497 -1.23E-06 
Santa Rosa 1 0.9522 0.0952 0.9522 0.1759 -1.12E-06 
Santa Rosa 2 0.9590 0.0767 0.9590 0.1491 1.84E-06 
EN A 
Ruiz Cortines 0.9999 2.67E-006 0.9999 0.0009 7.76E-07 
Batequis 0.9999 1.80E-008 0.9999 6.36E-05 7.25E-07 
Santa Rosa 1 0.9999 -4.01E-009 0.9999 7.49E-05 -8.33E-07 
Santa Rosa 2 0.9999 -2.87E-008 0.9999 8.45E-05 -7.88E-07 


R*: coeficiente de determinación; RMSE: raíz cuadrada del cuadrado medio del error (kPa); MBE: error medio sesgado; m: pendiente; b: intercepto; 
T,,, y: temperatura mínima; T”,, : temperatura máxima; AR: humedad relativa mínima; HR. : humedad relativa máxima: E; escenario ¡-ésimo. 


n máx máx 


Cuadro 3. Estadísticos de prueba de ajuste entre los datos observados de presión real de vapor de agua y los estimados 
por los modelos convencionales para las diferentes estaciones. 


Mejores modelos convencionales m b RS RMSE MBE 
Modelo 3 
Ruiz Cortines 0.9239 -0.0857 0.9349 0.3088 -0,2352 
Batequis 0.9307 0.0030 0.9153 0.2600 -0.1255 
Santa Rosa 1 0.9405 -0.1203 0.9203 0.3300 -0.2388 
Santa Rosa 2 0.9816 -0.0532 0.9124 0.2408 -0.0876 
Modelo 5 
Ruiz Cortines 1.0325 -0.0664 0.9349 0.2150 -0.0024 
Batequis 1.0402 0.0324 0.9153 0.2712 0.1069 
Santa Rosa 1 1.0508 -0.1048 0.9204 0.2521 -0.0036 
Santa Rosa 2 1.0966 -0.0293 0.9124 0.3007 0.1511 
Modelo 6 
Ruiz Cortines 0.9748 0.2772 0.9769 0.2569 0.2276 
Batequis 0.9551 0.3219 0.9793 0.2646 0.2386 
da Santa Rosa 1 0.9738 0.2670 0.9751 0.2494 0.2149 
3 Santa Rosa 2 0.9791 0.3262 0.9663 0.3176 0.2873 
a Modelo 7 
2 Ruiz Cortines 0.9827 0.2762 0.9755 0.2714 0.2422 
El Batequis 0.9640 0.3198 0.9785 0.2781 0.2531 
E Santa Rosa 1 0.9816 0.2659 0.9739 0.2636 0.2293 
e Santa Rosa 2 0.9904 0.3215 0.9653 0.3336 0.3036 
E Modelo 8 
be Ruiz Cortines 0.9786 0.2689 0.9767 0.2563 0.2267 
dE Batequis 0.9590 0.3136 0.9791 0.2637 0.2377 
z Santa Rosa 1 0.9775 0.2588 0.9752 0.2489 0.2141 
3 Santa Rosa 2 0.9839 0.3171 0.9659 0.3178 0.2869 
3 R?: coeficiente de determinación; RMSE: raíz cuadrada del cuadrado medio del error (kPa); MBE: error medio sesgado; 11: pendiente; b: 
E intercepto. 
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modelos 6, 7 y 8 presentan R? muy cercanos, y 
estos últimos tienen valores de R? más altos (en 
unas centésimas) que los modelos 3 y 5. 

Comparando los RMSE del modelo 3 (T... 
- 2) contra los modelos del 5 al 8, el modelo 
3 es mejor (por unas centésimas) en todos los 
casos en Batequis y Santa Rosa 2; además, al 
contrastar el modelo 7 con los modelos 5, 6 y 8, 
estos últimos resultan ser mejores en las cuatro 
estaciones analizadas. Generalizando, de los 
modelos convencionales, el mejor modelo es 
el 5, al considerar los estadísticos de ajuste; 
le seguirían los modelos 3, 6 y 8, y por último 
el modelo 7, pero en la práctica todos ofrecen 
una buena alternativa y su selección está en 
función de la información disponible, ya que 
las variaciones que presentan en sus RMSE 
entre ellos son mínimas. 

Se puede decir que los estadísticos RMSE 
y MBE de los dos modelos de RNA tienden a 


estar más cerca de cero que los RMSE y los MBE 
de los modelos convencionales; en particular, 
el modelo 5 de los convencionales muestra ser 
también un buen modelo de ajuste de la presión 
real de vapor. Sin embargo, las R? para modelos 
convencionales y modelos de RNA mostraron 
tener valores más o menos similares para todos 
los modelos, excepto para los modelos de RNA 
feedforward backpropagation y radial basis function 
con el mayor número de variables de entrada, 
que obtuvieron las R? más altas. 


Validación de los modelos 


En las figuras de la 6 a la 8 se observa que las 
diferencias entre los valores estimados con 
el modelo 1 E4 y los observados en los 229 
días están alrededor del valor cero (figura 6). 
El que presenta las diferencias más grandes 
es el modelo 3 (figura 7); el modelo 5 (figura 
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Figura 6. Diferencias del modelo 1 con los 229 días de datos en las diferentes estaciones. 
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Figura 7. Diferencias del modelo 3 con los 229 días de datos en las diferentes estaciones. 


8) presenta una condición intermedia con 
respecto a los dos anteriores. 

Al observar los RMSE y MBE del cuadro 4, 
se define que el modelo 2, en la mayoría de los 
escenarios y estaciones, es mejor que el modelo 
1, son la excepción E2 en la estación Batequis 
y el El en la estación Santa Rosa 2. Pero por 
las diferencias tan pequeñas que se presentan 
(décimas o milésimas) entre los dos modelos 
en todos los escenarios tanto en el ajuste como 
en la validación, se pueden utilizar de manera 
indistinta. 

De los modelos restantes, el modelo 3 es 
el que mayor error presentó en la validación, 
y el que tiene el menor es el modelo 5, que 
mostró consistencia para los dos conjuntos de 
datos usados (ajuste y validación); le siguen los 
modelos 6, 7 y 8 y estos fueron mejores en Santa 
Rosa 1. 


El modelo 2 E4, en el ajuste y en la 
validación, fue el que mejor ajustó a los datos 
observados; sin embargo, los entrenamientos 
que contaron con una sola variable de entrada 
también presentaron un buen ajuste de los 
datos observados, lo cual indica simplicidad y 
mayor interés práctico. No hay que descartar 
los modelos empíricos convencionales, como 
el modelo 5, ya que los valores RMSE son 
menores que los de los modelos 1 y 2 en 
sus escenarios El y E2. Fue la excepción la 
estación Santa Rosa 1, en este último escenario 
(E2) en ambos modelos (1 y 2). Además, 
el conocimiento de dichos modelos en sus 
variables describe al fenómeno en cuestión 
físicamente, y se refleja en los parámetros que 
integran a la ecuación, esto sirve de apoyo 
para conocer las variables de entrada en el 
entrenamiento de las RNA para estimar la 
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Figura 8. Diferencias del modelo 5 con los 229 días de datos en las diferentes estaciones. 


Cuadro 4. RMSE y MBE obtenidos en la validación, de las diferentes estaciones y con los diferentes modelos. 


Estaciones 
Modelos Ruiz Cortines Batequis Santa Rosa 1 Santa Rosa 2 
RMSE MBE RMSE MBE RMSE MBE RMSE MBE 
1, El 0.2426 -0.1480 0.2542 -0.1247 0.2716 -0.1569 0.3557 -0.2712 
EZ 0.2291 -0.1567 0.2487 =0/1855 0.2553 -0.1618 0.3924 =0:3127 
1, E3 0.2139 -0.0956 0.1728 -0.0402 0.2307 -0.1169 0.2629 -0.1705 
1 E4 0.0212 -0.0034 0.0253 -0.0023 0.0195 -0.0003 0.0972 -0.0481 
2, El 0.2411 -0.1460 0.2469 -0.1117 0.2656 -0.1470 0.3713 -0.2780 2 
a 
2, E2 0.2287 -0.1552 0.2519 -0.1293 0.2464 -0.1525 0.3705 -0.2847 2 
2, E3 0.2105 -0.0950 0.1655 -0.0237 0.2264 -0.0925 0.2327 -0.1361 E 
2, E4 0.0011 -0.0001 0.0001 -1.11B-05 0.0002 -2.76E-05 0.0003 -0.0001 = 
3 0.4307 -0.3824 0.3173 -0.2315 0.4484 -0.3845 0.3873 -0.3217 he 
5 0.2291 -0.1030 0.2301 0.0532 0.2707 -0.0984 0.2143 -0.0189 E 
6 0.2357 0.2013 0.2549 0.2297 0.2120 0,1794 0.2861 0.2556 le 
Y 0.2542 0.2187 0.2747 0.2493 0.2306 0.1968 0.3089 0.2778 E 
8 0.2372 0.2022 0.2563 0.2309 0.2137 0,1805 0.2891 0.2580 3 
Ei: escenario ¡-ésimo; RMSE: raíz cuadrada del cuadrado medio del error (kPa); MBE: error medio sesgado. 3 
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presión real de vapor de agua. En el ajuste, 
el modelo 3 aproximó muy bien a los datos 
observados de las estaciones Batequis y Santa 
Rosa 2, pero en la etapa de validación presentó 
los valores más grandes de RMSE y MBE en 
las cuatro estaciones. 

En el ajuste, el modelo 3 fue mejor que el 
modelo 5 en Batequis y Santa Rosa 2 debido a 
que el promedio de la humedad relativa para el 
total de datos que se utilizaron (1 484) en cada 
estación fue de 66.4 y 65.7%, respectivamente. 
En las otras dos estaciones (Ruiz Cortines y 
Santa Rosa 1) donde fue mejor el modelo 5, 
la humedad relativa fue de 70.3 y 70%. En la 
validación, el modelo 5 fue mejor que el 3 en 
todas las estaciones, ya que la humedad relativa 
de 229 datos que se utilizaron en Ruiz Cortines, 
Batequis, Santa Rosa 1 y Santa Rosa 2 fueron 
de 74.2, 70.4, 74.7 y 71.7%, respectivamente. 
Doorenbos y Pruitt (1980) indican que una 
región con humedad relativa promedio igual 
o mayor al 70% se clasifica como una zona 
con humedad alta; en el proceso de ajuste, 
dos estaciones cumplen con esta condición y 
en la validación, las cuatro. Por tal motivo, el 
modelo 5 fue mejor, lo que concuerda con lo 
encontrado por Allen et al. (1998) y Cai et al. 
(2007), que en zonas húmedas o en localidades 
donde el cultivo está bien regado, la T,. =T...... 

Con base en los resultados obtenidos, los 
modelos 1 y 2, en sus escenarios E3 y E4, se 
pueden utilizar en trabajos de investigación 
con modelos de simulación de cultivos para 
pronóstico de la evapotranspiración en tiempo 
real en sistemas de riego de alta frecuencia. Los 
otros escenarios de los modelos 1 y 2, el modelo 
5 y los modelos 6, 7 y 8, se sugiere se utilicen en 
la planeación y operación del riego de cultivos 
donde se requiera cada tercer día o una vez a la 
semana. El modelo 5 demostró en condiciones 
de alta humedad (zona húmeda o regiones con 
cultivo bien regados) que es buen estimador de 
la presión de vapor de agua real. En regiones 
secas o de baja humedad relativa, se necesitan 
calibrar los modelos 3 y 4, para definir cuál es 
el valor que se le va a restar a la temperatura de 
punto de rocío, para que estime bien la presión 
de vapor de agua real. 


Conclusiones 


De los modelos convencionales, el que mejor 
aproximó a los datos observados de manera 
global fue el modelo 5 (T,. ). Delos dos modelos 
de RNA (feedforward backpropagation y radial basis 
function), el que realizó mejor aproximación a 
los datos observados de presión real de vapor 
de agua es la RNA radial basis function E4 
(con las variables de entrada de temperatura 
máxima y mínima, humedad relativa máxima y 
mínima). De manera general, se puede concluir 
que las redes neuronales artificiales (RNA) 
muestran ser unas herramientas robustas en la 
estimación de la presión real de vapor de agua, 
en sus escenarios E4. 
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Abstract 


CERVANTES-OSORNIO, R., ARTEAGA-RAMÍREZ, R., VÁZQUEZ-PEÑA, M.A., 
OJEDA-BUSTAMANTE, W. € QUEVEDO-NOLASCO, A. Comparison Models for 
calculating actual vapor pressure. Water Technology and Sciences (in Spanish). Vol. IV, 
No. 2, April-June, 2013, pp. 37-54. 


Although real vapor pressure is a fundamental variable in the calculation of the 
evapotranspiration of crops, a component of the hydrological cycle, it is difficult and costly 
to measure directly. Consequently, in practice calculations are based on temperature data 
and psychometric relationships. The objective of the present work was to compare different 
conventional methods used to calculate real vapor pressure with calculations performed using 
two types of artificial neural networks-Feedforward backpropagation and radial basis function. 
Meteorological data were used from four weather stations in Irrigation District 075, located in 
Valle del Fuerte, in northern Sinaloa State, Mexico. The results indicate that the radial basis 
function type of artificial neural network (Scenario E4) was the best method to calculate real 
vapor pressure. 


Keywords: atmospheric humidity, vapor pressure deficit, artificial neural network. 
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